Transient detections and other real-time data processing from 
MASTER- VWF wide-field cameras. 



Evgeny Gorbovskoy 1 , Kirill Ivanov 4 , Vladimir Lipunov 1 , Victor Kornilov 1 , Alexander Belinski 1 , Nikolaj 
Shatskij 1 , Dmitry Kuvshinov 1 , Nataly Tyurina 1 , Pavel Balanutsa 1 , Vadim Chazov 1 , Artem Kuznetsov 1 , 
Petr Kortunov 1 , Andrey Tlatov 2 , Alexander Parkhomenko 2 , Vadim Krushinsky 3 , Ivan Zalozhnyh 3 , 
Alexander Popov 3 , Taisia Kopytova 3 , Sergey Yazev 4 , Alexander Krylov 1 

1 Moscow State University, Sternberg Astronomical Institute, 119991, 13, Univeristetskij pr-t, Moscow, Russia 

2 Kislovodsk Solar Station, 357700 p.o. Box 145, 100, Gagarina st., Russia 

3 Ural State University, 620083, 51, Lenina pr-t, Ekaterinburg, Russia 

4 Irkutsk State University, 664003, 1, Karl Marks st., Irkutsk, Russia 

Construction of robotic observatories has developed into an important and thriving fields of astronomy. 
Their large field of view combined with the capability to be pointed at any direction make robotic astro- 
nomical systems indispensable for tasks involving searches for transients like grb, super novae explosions, 
novae etc., where both the time and direction of the search are impossible to predict. This paper de- 
scribes prompt GRB observations made with MASTER- VWF wide-field cameras and the methods of 
image analysis and classification of transients used for real-time data processing. During seven months of 
operation six synchronous observations of gamma-ray bursts have been made using MASTER VWF facil- 
ities deployed in Kislovodsk and Irkutsk. In all cases high upper limits have been obtained (see Table [ 
and Fig. [Q). 



1 Introduction 

1.1 Very wide field telescopes. 

Currently, several dozen robotic telescopes are operated worldwide, each of them usually meant for certain 
dedicated tasks. A special case of these instruments are super-wide field telescopes with fields of view 
ranging from hundreds to thousands square degrees. It is impossible to make a very wide field system 
with a large aperture using now available CCDs and optical systems, and that is why the apertures of 
very wide field cameras are limited to 5-15 cm. All very wide field cameras are designed to record various 
optical transients. Most of the optical transients are of circumterraneous origin — meteors and satellites, 
however, very wide-field cameras may also provide valuable data on the prompt emission of cosmological 
grbs. A striking example is the incredibly successful observation of GRB080319B that was performed 
within the framework of the Russian-Italian experiment TORTORA [5j and Polish experiment " Pi of the 
Sky" [11]. All transients are highly variable and therefore snapshots must be taken at intervals of several 
seconds or even less. Hence solutions capable of operating without dead time (such as the time needed for 
ccd readout), like MASTER VWF pQ, TORTORA and others (see table [T]), have an important advantage 
over other systems. 

The operation of any automatic astronomical facility is supported by a dedicated software package, 
which also performs primary image processing and data storage. Which software to use for an automatic 
astronomical system depends on the particular task the facility is meant to perform, and multi-purpose 
solutions are therefore unsuitable in such cases. Classical universal techniques are even less suited for 
super wide field observations, because of the problems due to the field curvature and very large data 
amount. 



1.2 MASTER - Mobile Astronomical System of the TElescope Robot 



MASTER was the first robotic telescope of this kind constructed in Russia. It was built and tested by 
students, post-graduate students and lecturers of SAI MSU in 2002 (http:/ /observ. pereplet.ru). MASTER 
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Table 1: Comparative characteristics of super- wide field systems 
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85x70 
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Fig. 1. Summary results of synchronous observations of a gamma-ray burst. Dependence plotted as a function 
of gamma-ray fluence in the soft (15 — 150keV) and hard (8 — 1000/ceV) energy interval. 



has evolved into a broad network of telescopes spread across the entire country. A detailed description of 
the MASTER facility can be found in [2] . Below we discuss the basic properties of MASTER wield field 
cameras (M ASTER- VWF). 

The main task of the MASTER- VWF experiment is to perform continuous all-sky monitoring with 
the aim to detect all objects absent in the available astronomical catalogues. This includes, in particular: 

• Prompt detection of grb emission synchronously with gamma-ray space observatories. 

• Search for star-like transients of unknown origin and orphan bursts. 

• Detection of meteors and determination of their basic parameters: brightness (photometry), velocity 
(astrometry) , and altitude of combustion in the atmosphere (triangulation). 

• Detection of satellites and space debris and determination of their basic parameters - astrometry, 
velocity, brightness, and altitude above the earth's surface. 

For these purposes, a completely robotic MASTER VWF4 facility with four wide-field cameras and a 
total field of view of 4000 square degrees has been deployed at the MSU Caucasian Mountain Astronomical 
Observatory near Kislovodsk(2075m). Another facility — MASTER VWF2 (two cameras and 2000 square 
degrees field of view) — has been deployed near Irkutsk. In Kislovodsk cameras are set on parallactic 
mounts in pairs 702 meters apart (see fig. [2]). Each mount has an automated dome and two fast and 
powerful (11 Megapixels) Prosilica GE 4000 CCD cameras with Nikkor 50mm (//1.2) lenses. The camera 
is capable of continuous imaging with a speed of one frame in 0.2 to 60 seconds with no time gaps. The 
702-m baseline allows the facility in Kislovodsk to determine the parallaxes of circumterraneous objects 
and reconstruct their altitudes from the results of astrometric reduction. 

All telescopes operate in a fully autonomous mode: at night, if the weather conditions permit (a 
"Boltwood Cloud Sensor" provides the current weather report updated every three seconds), the cameras 
begin to automatically survey the sky. The frames taken are then processed and stored in a special data 
storage device. Because of the very large data flow (up to 700 Gb per night), the images can be kept 
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Fig. 2. The equipment of the MASTER facility deployed at the MSU Caucasian Mountain Astronomical Obser- 
vatory near Kislovodsk. From the left to right: (1) the East and West cameras of the MASTER- VWF4 northern 
installation; (2) the main tower and shelter of the southern MASTER- VWF4 installation, and (3) the main tower 
of the MASTER- 2 telescope [2] 



only for several days and then most of them have to be erased, except for some images with objects of 
special interest like meteors, candidate optical transients, and satellites. We also keep the images taken 
simultaneously with grbs. The survey is halted and the telescope dome is closed every time when the 
weather deteriorates or at dawn. The operation of the entire MASTER VWF4 facility is controlled by 
seven computers (see Fig. [3]): four computers to control the operation of the Prosilica GE 4000 CCD 
camera and perform source extraction; two computers to control the operation of the dome and mount; 
one big quad-core server to process and store the data obtained and publish it in the Web. 

The facility is a complex engineering system, which needs sophisticated and powerful software to con- 
trol it, perform observations, and data processing. The data flow can be as high as several Terabytes/day 
(in the case of very short exposures on the order of 0.2s), requiring highly optimized software, powerful 
data storage facilities, and interfaces. With 5-second long exposures the limiting is 11.5-12 on each frame 
depending on weather conditions. Thus every image contains about 10-15 thousand objects, resulting in 
a data flow that is impossible to store and has to bee processed in real time. 

A special software package has been developed for this purpose. It allows the coordinates of all 
objects and their photometric characteristics (profile, magnitude) to be determined in real time, and is 
also capable of finding and analysing optical transients. MASTER VWF can measure the coordinates of 
the objects with an error smaller than 10 arcsec. The system is supplied with a daily updated database 
of the ephemerides of artificial satellites and a program developed to identify the satellites in real time. 
The program is capable of finding new artificial satellites, to observe the destruction of known artificial 
satellites and not only to control space debris, but also observe the process of its formation. Sometimes 
unknown or lost objects are recorded. For example, on October, 20th, 2008 MASTER VWF recorded the 
debris of the American military satellite USA114Deb (see fig. 

A total of six synchronous observations of gamma-ray bursts have been made by MASTER VWF in 
Kislovodsk and Irkutsk during seven months of their operation. In all cases high upper limits have been 
determined (see tabl 14.31 and fig. [1]). 

In the following chapters we describe these and other results in more detail. 

2 Primary reduction in wide fields: astrometry and photometry 

Let us consider the methods of astromeric and photometric reduction of a wide field in more detail. The 
identification of the field of view and calculation of the coordinate grid for the image is based on the 
table of Cartesian image-plane coordinates (X,Y) produced by program used for the extraction of objects 
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Fig. 3. Schematic diagram of the layout of the equipment in the MASTER VWF4 experiment. 



(SExtractor) and the the Tycho-2 astrometric catalogue, which contains exact coordinates for most of 
the stars down to a limiting magnitude of ll m .5. To perform this task, the program should meet high 
requirements in terms of speed of processing (10-15 seconds and less) for a field containing about 10000 
stars; ensure high probability of identification (more 99 %); be applicable to various fields ranging in size 
from from narrow (tens of square minutes) to super- wide (thousand square degrees) fields; allow for a 
certain freedom in the telescope parameters, because the focal length and other parameters may change 
in the course of the year due to the variation of the outside temperature, and be capable of operating 
with approximate initial midfield coordinates provided that they are within the field of view. 

All processing in MASTER experiments is performed under Linux OS and hence no ready-made 
software was likely to be found. The current implementation of the "astrom" program is based on the 
algorithm of pairing similar triangles of bright stars in the image and in the catalogue. This program 
works perfectly with the most different images from those taken under poor weather conditions and 
containing a few dozen stars (typical situation for alert observations) to images in the Milky- Way fields 
taken during a good observing night and containing up to 100000 objects. The root-mean-square error 
of the measured coordinates ranges from 1 / \ of arcsecond (scale 2.1" /pixel) for most of the images taken 
with the basic telescope of the MASTER system to about 7-8 arcseconds for very wide-field cameras 
(scale 36" /pixel) (see fig. |2j). 

The system may also operate in a "fast" astrometry mode to process consecutive frames of the same 
area taken with wide-field cameras. In this mode the software uses a priori information about the field of 
view based on the previous frames, allowing the operating time to be reduced to less than 0.2 seconds. 
There is also a dedicated mode for more accurate determination of the coordinates of special objects 
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Fig. 4. Debris of the American military satellite USA114Deb found via the web interface of the MASTER 
system. The top and bottom rows of the images are taken by the Northern and Southern cameras, respectively. 
The parallactic shift is apparent event to an unaided eye. The measured parallax implies a height of 4000 =b 200 
kilometers 



(gamma-ray bursts, supernovae, asteroids, meteors, and others) with the local coordinate adjustment 
computed only in the field of the object searched using 10-25 nearest stars with well accurate coordinates. 
Our experience shows that this mode provides a factor of two improvement of the accuracy of the inferred 
coordinates. 

Our program for extracting objects from images is based on the SExtractor software package written 
for the TERAPIX [9j project. SExtractor is written in the C language and offers sufficiently high per- 
formance. We adapted this program to the task of reducing images taken with MASTER, and describe 
our version of the program in the following chapter. The main parameter used by the extraction program 
is the threshold signal-to-noise ratio. This value directly affects both the program execution time and 
the accuracy of the results obtained. We adopt S/N=2 for primary selection, whereas at the subsequent 
stages of processing we select from the initial list only the objects with higher S/N values. 

2.1 Photometric calibration 

Practically all wide-field systems operate without filters. MASTER VWF is not an exception. This poses 
the problem of the calibration of the photometric band of the measurements. It is clear that our instru- 
mental photometric band is limited to a certain extent by the CCD sensitivity curve and lens transmission 
curve. Our task is therefore to determine instrumental magnitudes based on the magnitudes of the refer- 
ence stars adopted from photometric catalogues. 

To this end, colour-colour diagrams are analysed, which relate the difference between the instrumental 
magnitude and the magnitude in any of the modelled bands (or combinations thereof) to a chosen colour 
index for each detected star. The closer is the modelled band to the instrumental band, the more horizontal 
appears the the cloud of points on the diagram. Fig. [6] shows the corresponding colour-colour diagrams 
for the MASTER VWF-4 system consisting (in its photometric part) of a Nikkor 50 mm f/1.4 lens 
and Prosilica GE4000 CCD. The instrumental photometric band can be described fairly well by the 
Vtycho2 filter. The absolute accuracy of photometric measurements made with wide-field cameras is, 
on the average, of about Am a bs ~ 0.25 — 0.35 m . Note that the relative accuracy of measurements is 
substantially better with errors as low as Am re i 3 — 5%. 

2.2 Automatic determination of the limiting magnitude of an image 

Automatic determination of the limiting magnitude for the images taken is an important part of the 
operation of modern automatic sky surveys. An automatic survey must naturally include automated 
quality control. The automatically determined limiting magnitude may serve as some kind of objective 
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Fig. 5. Average error (the upper curve) and the number of identified stars (the bottom curve) as functions 
of the degree of the polynomial used. The dependences are shown for three different optical systems. Red (left) 
Rigter-Slefogt (MASTER) with a field of view of 6 n °, Nikkor 50 mm f/1.2 (middle) with a field of view of 1000 n °, 
and Nikkor 50 mm f/1.4 (right) with a field of view of 1000n°- The curves are based on arbitrary frames from 
various sky areas. A third to fourth degree polynomial approximation can be seen to provide optimum solution 
(95-99 percent of stars are identified with the minimum average error). A more detailed analysis shows that 
the third-degree fit yields slightly better results. That is logical, because most of the known aberrations (e.g., 
distortion) depend on the third power of radial distance. 



quality assessment for the images obtained. Furthermore, the determination of the limiting magnitude is 
also needed for gamma-ray burst observations. 

The task of the determination of the limiting magnitude can be addressed formally by setting the 
limiting magnitude of the frame equal to the magnitude of the faintest object. This approach, however, 
has a number of shortcomings. First, when the system is operated with a low extraction threshold most of 
the faintest objects are noise artifacts. If the sky is covered with clouds, no real faint objects can be seen 
and practically all objects at the visibility limit are chance fluctuations of cloud brightness. Moreover, 
because of the large number of such false objects some of them may be identified with real stars and 
hence the situation cannot be remedied by restricting the list to objects identified with catalogued stars. 
Second, this approach allows one to estimate the limiting magnitude only within a local area and the 
quantity inferred does not characterise the entire image and is therefore unsuitable for quality control. 

These problems can be addressed as follows. All stars of the reference catalogue that are located 
within the field considered are subdivided into magnitude bins of fixed width. For each bin the fraction 
of identified stars to the total number of catalogued stars is computed, and the corresponding histogram 
is drawn (see fig. [7j). One can see that the ratio decreases sharply when the average magnitude of the bin 
approaches the limiting magnitude. In practice, it is convenient to adopt a certain threshold ratio and 
define the limiting magnitude as the magnitude at which the curve reaches this threshold. The thresholds 
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Fig. 6. colour-colour diagrams are used to calibrate instrumental magnitudes. On these plots m — 2.5lg(Fi ns t) 
and the zero point is of no importance and cannot be determined at the this stage. Note that the width of the 
sequences is real and due to the fact that stars have different spectral types, and not due to measurement errors. 



of Gum = 25% and au m = 50% have been shown empirically to correspond to objects with 5 /at~9 — 11a 
and s / n ~ 14 — 15a, respectively. 

Note that for robust operation of the algorithm the limiting magnitude of the reference catalogue 
should be fainter than that of the frame studied. Otherwise the method becomes inapplicable because 
the catalogue lacks lack stars of the given magnitude. However, the USNO-B1 catalogue is sufficient for 
most of the modern wide-field surveys, whereas very wide-field surveys can be reduced successfully even 
with the TYCHO-2 catalogue (and this is very convenient). 



3 Automatic classification and analysis of astronomical transients. 
3.1 Extraction of objects and data transfer to the server. 

The search for and classification of transients starts once the primary reduction of the image is finished. 
All transient events in frames can be subdivided into two classes: starlike objects and strips, which differ 
substantially in the methods used for primary reduction. All rapidly moving objects are usually classified 
as strips: low-altitude satellites and meteors. All transients of stellar or extragalactic nature (novas and 
supernovas, orphan bursts, stellar flares) and geostationary satellites (see, for example, figIT8]) are starlike. 
Starlike objects are much easier to process, because the program of object extraction and instrumental 
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10 12 14 16 18 20 10 12 14 16 18 20 

Magnitude Magnitude 

Fig. 7. The number of identified stars as a function of magnitude. The straight lines indicate the an m = 50% 
and Gum — 25% thresholds. Averaging was made within 0.5 m and 0.1 m bins in the right and left histograms, 
respectively. The right (coarser binned) histogram is smoother, resulting in reduced probability of accidentally 
reaching the threshold at a bright magnitude. It is therefore convenient to determine the limit on the coarsely 
binned histogram and then refine it on the finer-binned curve. 

photometry (SExtractor, described above) has been developed to analyse such objects. Therefore all main 
parameters of star like objects (flux and magnitude, FWHM, Cartesian frame coordinates (X, Y), and 
semiaxes (a, b)) are determined very reliably. 

Extraction of strip-shaped objects is fraught with serious problems. First, a fine adjustment of S Ex- 
tractor's additional internal parameters is necessary. Second, a change of these parameters strongly affects 
the image processing time. Furthermore, strip-shaped objects contribute appreciably to the background 
within the small area suitable for identification of stars. Fine adjustment of all parameters allowed us 
nevertheless not only to identify strips in the image, but also achieve high speed of operation (half expo- 
sure time) (see fig IE). However, SExtractor algorithm cannot identify strips in the frame, to say nothing 
about performing their high-quality photometry. To this end, a special program had to be developed, 
which we describe below. 

To reduce the server load, identification of objects in the image can be physically performed not only 
by the server, but also by the control computer of each particular camera (see fig. [3). Note that the 
resulting catalogue of detected objects should have higher priority for the transmission to the server 
compared to the images that are to be stored in the database. 

3.2 Primary reduction of astronomical transients. 

The next stage after identification of objects is their astrometric and photometric reduction, which we 
describe in detail in the previous section. We thus obtain, after identification, astrometric and photomet- 
ric calibration, a catalogue of objects identified in the frame together with calibrated coordinates and 
photometry. All unidentified objects (identification is performed using both coordinate and photometric 
criteria: usually AS — 20 — 30" and Am = 1.5 m ) are viewed as possible optical transients. This approach 
allows us to identify stellar flares as transients and to a certain extent prevent spurious identifications 
with fainter and hence more numerous stars. To filter noise, objects that could not be identified with the 
catalogue are correlated with the objects of several (usually 2-3) previous sets. An object is considered to 
be a candidate transient if it could not be found either in the catalogue or in in the previous frames. All 
candidate transients are then passed through a number of special filters. For example, if the semimajor 
axis of the object, which should typically be of about b ~ 0.7 — 2 pix, exceeds 6^4 — 5 pixels, it must 
be a trail of a cloud with a probability of 95 %. If fwhm < 1 pixel, we must be dealing with a hot pixel 
or a cosmic-ray hit. We also set for all transients the trigger threshold, which is always higher than the 
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Fig. 8. A 5 o x 9o image with a bolide. The black squares indicate objects extracted from the frame. The side 
of the square corresponds to the fwhm of the object. The left-hand figure shows a frame with objects extracted 
using standard SExtractor parameters (with the extraction threshold of s /n > 1.5) with no filters applied. Ax 
a result, the program identified all starlike objects, heaps of faint garbage, hot pixels, and simply random noise, 
i.e., virtually everything except the bolide. The right-hand frame shows the result of extraction made using the 
currently adopted refined parameters. Several bright elongated objects are extracted on the bolide image and the 
noise level is appreciably lower. 

extraction threshold. Currently, the trigger threshold is s /at > 10a. We consider a candidate found on 
two or more consecutive frames to be a boda fide transient. 

Preliminary classification of transients is performed after the objects are passed through all filters. 
First, all objects are searched for possible identifications in the catalogue of artificial satellites and all 
identified objects are included into the database of satellites. The remaining ones are subdivided into 
starlike and strip-shaped objects. The object is classified as a strip-shaped if a/b > 3, otherwise it is 
considered to be star like (here a and b major and minor semiaxis). All starlike objects are then passed 
through the coincidence circuit, and strip-shaped objects are. Strips are provisionally subdivided into 
meteor candidates and unknown artificial satellites. If each of the three successive frames contains a 
candidate transient such that the three objects lie on a straight line, all of the them are considered to be 
the images of an artificial satellite. The remaining objects are considered to be meteor candidates. 

3.3 Coincidence scheme of and determination of heights. 

All the procedures described above are performed for each of the four cameras in real time mode. Real time 
processing ends here if the cameras operate in uncoupled (alert) mode. In the case of survey observations, 
when two cameras point to the same sky area, all starlike transients are passed through the coincidence 
scheme, and all strip-shaped objects have their heights determined (see fig. I§1 [TP]). To improve the accuracy 
of identification, special transformation is applied to the frame taken with the Northern camera to convert 
it to the coordinate frame of the Southern camera. This transformation is based on the same principles 
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Fig. 9. MASTER VWF4 frames of a 5° x 5° sky area taken with the Northern and Southern cameras. The 
parallax of the meteor is immediately apparent. Triangulation yields a height of H meteor = 72 ± of 2 km. 
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Fig. 10. Parallax of an artificial satellite. Here we show the 4° x 4° frames of a sky area taken with the MASTER 
VWF4 Northern and Southern cameras. Triangulation yields a height of H sate ii ~ 4500 km. Given the accuracy 
of astrometric measurements for strip-shaped objects, cr^^from ~ 20", and the baselength of Al = 702 m, the 
maximum height that can be determined via triangulation is H ma x = Al * arcsin{a l a™trorn) ~ 10000m. 



as astrometric reduction. It takes little time to compute (~ 0.2 seconds) and improves the accuracy (and 
hence reduces the correlation radius and the number of chance identifications) by more than a factor of 
2.5 compared to the accuracy of the results obtained in the case of identification by celestial coordinates. 
Should a real starlike transient be found, an alert can be immediately issued to larger (45-cm) MASTER 
telescopes (located near Moscow, in Kislovodsk and in the Urals) to study the event in more detail. 
However, this mode is currently in the stage of development, mostly because of the debris of high-orbital 
satellites, which often contaminate the frames (see fig. IH IT8|) . 
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3.4 Analysis of strips. 

As we pointed out above, SExtr actor computes only the approximate position of a strip, breaking it into 
several bright elongated objects. To perform complete photometry of a strip, determine the coordinates 
of its centre and ends and other parameters, a special program is used dedicated to the analysis of strips 
on astronomical images. To prevent server overload, this program is run only when cameras are idle (in 
the afternoon or during bad weather). The program starts with the coordinates of any location inside the 
strip determined at the previous stage. The working area of the given size (usually 512 on 512 pix, which 
is sufficient for the analyses of any satellite or meteor trails) is chosen around this initial point. The size 
of the working area is automatically increased when a very bright and long bolide is to be analysed. 

Determination of background. 

To determine the background, we subdivide the working area into subareas of the given size. For each 
subarea the median and standard deviation are computed twice: the first time over almost all pixels in the 
given subarea, and the second time more precisely, after discarding the pixels that deviate strongly (by 
more than 6- 10<j) from the preliminary background estimate. We number all these subareas — there are 
a total of N zone of them (currently N zone = 256) — and compute the background and standard deviation 
in each subarea as 



where Bg (i,j) is the background at the point with coordinates (i,j); Med (k) is the median in k-th area 
used to compute the background; £ is the standard deviation within the corresponding subaera, and r& 
is the distance between the point (i.j) and the centre of k-th area. Figures [TO 3-4 show the corresponding 
background maps. 

Extraction of strips. 

To extract a strip, we first compute the map of bright pixels exceeding the background level by n a • j) 
or more. The optimum threshold has been shown empirically to be n a = 1.5. The corresponding points 
are shown by black shade on the map (see fig [TTJ2). 

A recursive algorithm of the search for connected area is then used starting at the centre found by 
SExtractor. Once the connected area is found, the main parameters of the object (the position of the 
centre, a, b, and the angle between the axis of the object and the X axis) are determined. The standard 
deviation au ne from the straight line drawn at angle 6 through the centre of the object is computed if 
the object found can be classified as a strip (a/b > 3). The recursive algorithm of searching for connected 
area is then started again, this time within 2-3<r^ ne of the straight line just mentioned. The algorithm 
is allowed to jump over one to two faint pixels within 2au ne of the straight line. Thus the second pass 
yields a completely extracted strip including faint edges in the case of meteors. 

The program then attempts to find objects located on the extension of the strip line on the subsequent 
and previous frames. If such object can be found, the strip is classified as a satellite, otherwise, as a meteor. 
The brightness profiles along the strip are then computed both for satellites and meteors (see fig. [T2l ). 
Portions of images are kept for all transients. 

4 Synchronous observations of prompt grb emission. 

In this chapter we discuss synchronous grb observations. Note that during two years of the operation of 
the single-channel wide-field Kislovodsk camera we have made over 20 faster than 1 minute alert pointings 
published in the corresponding GCN-telegrams GRB 070224 (gen 6139), GRB 070223 (gen 6131), GRB 
070219 (gen 6113), GRB 061213 (gen 5915), GRB 061002 (gen 5677), GRB 060929 (gen 5657), and 
others [3]. 




(i) 




k 



k=l 
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Fig. 11. (1) Image of a meteor on a 3° x 3° frame (F(i,j)). 

(2) Extraction map. Black shade shows the points with F(i,j) > n a • white shade indicates the points 
included into the object under study. Grey shade indicates the background. 

(3) Background map - (Bg(i,j)) (bottom left) 

(4) Map of standard deviation (U(i,j)) (bottom right) 



4.1 The importance of synchronous observations. 

As is well known, optical grb glow can be roughly subdivided into two parts: prompt emission and 
afterglow. We do not discuss the afterglow phenomenon here, because it has already been thoroughly 
studied and it is not the target of wide field observations. Let us now consider prompt emission. 

Prompt optical grb emission is the emission that appears simultaneously with gamma rays detected 
by various space-based gamma observatories (Hete2, Swift, Fermi). Unlike afterglow, it bears information 
about the object that gives birth to grb, and not about its environment. There are two ways to detect 
prompt optical grb emission: synchronous observations and very quick alert pointing. To date, about 
dozen [8] optical prompt grb have been made and most of them were alert based. That is to say, when 
prompted by a quick trigger from a gamma-ray observatory, a robotic optical telescope points at the object 
within several tens of seconds after the corresponding gamma-ray event is detected. If the grb is long 
enough, its prompt emission is still detectable by this time. Another, less common type of observations, 
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are synchronous wide-field camera observations. Perhaps the most amazing example of the latter is 
the prompt GRB 0803 19B detection [5J. Synchronous observations are much more efficient than alert 
observations for at least two reasons. First, alert observations are subject to selection effects. The point 
is that an alert may result in the detection of prompt emission only if the grb lasts longer than 30-40 
seconds, because of the delay due to the time it takes to perform data processing onboard the space 
telescope, transmit the signal to the ground, and point the optical telescope. Hence alerts cannot be used 
to detect prompt optical emission from short grbs. It is significant that this kind of emission has never 
been observed and moreover, it has never had its upper limit measured. 

4.2 Energy calibrations of the MASTER- VWF4 photometric passband. 

To make any conclusions about the brightness of objects not just in terms of relative magnitudes, but also 
in terms of energy units, we must calibrate the M ASTER- VWF4 zero-point. Let us determine the energy 
flux (in the units if erg / cm 2 / s) produced by a zero-magnitude star (e.g., Vega). The M ASTER- VWF4 
band is close to the V band of the Tycho2 catalogue. Figure [13] shows the spectrum of Vega and the same 
spectrum folded with the lens passband and sensitivity curve of the CCD. Integration over the entire 
spectrum finally yields the flux produced within the MASTER- VW4 band by a zero-magnitude star 
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Fig. 14. Coadded images taken 30 seconds before (left), during (centre), and 30 seconds after (right) GRB081102. 
Each frame is a coadd of 12 images (six images from each camera) - 60 seconds exposure for each frame. Animation 
is available here: http://observ.pereplet.ru/images/GRB081102/grb_film.html 



F ve 9 a = / . u .du = 6A± 0.1 ■ 10~ 6 erg / cm 2 / s 



(2) 



4.3 Synchronous GRB observations with MASTER VWF very wide-field cameras 



MASTER VWF cameras perform nightly monitoring of the sky with the aim to detect optical tran- 
sients. These cameras may operate in two different modes. In the first mode, which is meant for searching 
for optical transients, two cameras are pointed to the same sky area in order to implement the coincidence 
scheme (see fig. [9j [T0|) . In the second mode the fields of view of the cameras are set apart as far as 
possible in order to provide a combined field of view of more than 4000 and 2000 deg 2 for the Kislovodsk 
and Irkutsk observatories, respectively. The second mode is dedicated for synchronous observations of 
gamma-ray bursts. 



The GRB081102 and GRB090424 gamma-ray bursts recorded by Swift spacecraft 

GRBs recorded by the Swift gamma-ray observatory satellite have very accurate coordinates, allowing 
the location of the burst origin on very wide-field cameras to be determined with the highest precision 
possible (i.e., to within less than one-pixel large error box). 

Long gamma-ray burst GRB081102 

(see fig. HH ) was recorded by Swift gamma-ray space observatory on November, 2nd, 2008 17:44:39.5 
UT [29]. Although it was a typical long (T90 = 63sec) gamma-ray burst, the telegram from Swift was 
issued only after 15 minutes. This delay prevented early observations by automated telescopes. However, 
MASTER VWF-4 in Kislovodsk covered the entire error box at the edge of its field of view. MASTER 
VWF-4 operated in the survey mode with aligned cameras, allowing the grb error-box to be recorded 
simultaneously in two channels. This area has been observed for two hours before and seven hours after 
the grb with no time gaps [T4j . 

Unfortunately, no optical counterpart could have been recorded from this burst (see fig. [I4j [T5|) . 
The lack of optical counterpart may have been due to extremely strong absorption toward the burst 
(jl]). However, a very high upper limit (for synchronous (!) observations with very wide-field cameras) of 
Vgr608ii02 < 13. m has been obtained, which is the highest synchronous upper limit in history [13]. 

Let us now compare GRB081102 with GRB080319B from which bright prompt optical emission V = 
5.3 m ([5]) has been recorded. Consider the F opt / 'F 1 ratio (where F opt and F 1 are the optical and gamma- 
ray flux, respectively) for both bursts. The gamma-ray spectrum is known for all bursts discussed in this 
paper and therefore we reduce F 1 to the flux in the 15 — 150/ceF energy interval for the sake of uniformity. 

Let us estimate F^ 081102 . The V-band extinction can be estimated by the following simple empirical 
formula [12] 
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Fig. 15. GRB081102. Hash (red) from the Northern and Southern cameras and from the two cameras simulta- 
neously inside the GRB081102 error-box and the signal from the V re f — 11. 5 m reference star (blue) 



Fig. 16. A 6x6 deg 2 field around the GRB090424 Swift XRT position. One-minute exposure images taken before, 
during and after Tqrb- Clouds can be seen on the images, resulting in the optical limit of V < 9.7 m on each 
frame. Animation available here: http://observ.pereplet.ru/images/GRB090424/grb-film.html 



r v = 5.2 • 10~ 22 N H (3) 

Given the known hydrogen column density toward GRB081102, Nh = 4.9 * 10 21 cm -2 [31J, we can 
use the inferred optical depth and Pogson's formula to estimate jp grb081102 = 2.3 • 10~ 6 erg/cm 2 [22] or 

rf 6081102 = 5.2 • 1(T 22 • 4.9 • 10 21 = 2.548 => 5m = 2.8 m (4) 



F opt° 81102 2 5i2 rnVega -( rn9rb081102 - 5rn ) • F Vega • T| r6081102 

1102 ^ 77,^6081102 (5) 



F°;r iW2 2.512°- 13 -°+ 2 - 8 • 6.4 • 10- 6 erg/s/cm 2 ■ 63s _ 1 
^6081102 < 2.3 • 10- 6 erg/cm 2 ~ 83 

The "less than" sign is used because only the upper limit could have been determined for grb081102. 
For GRB080319B the column density is Nh — 9.2 * 10 20 cm -2 [5j, i.e., five times lower than toward 
GRB081102 r f/ 60803195 = o.48 =; Sm = 0.5 m . Hence 



F grb080319B ^ 512 m Ve »°-(m* r6080319B -<5m) . pVega . T |r6080319B 

v ^r6080319B = ^6080319B (^) 



•7 ± 1 

-i#r6080319B 



F%r = 2.512°- 13 -°+ 2 - 8 • 6.4 ■ 10- 6 erg/s/cm 2 • 50s = 1 
F grb0803i9B ~ S2 5 • 10- 5 erg/cm 2 ~ 22 

It thus follows that the fraction of prompt optical emission for GRB081102 is at least four times 
smaller than for grb080319B. 
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Fig. 17. The most recent error box of GRB081130B a = 00^56 m 20 s 5 = +04°12 / , error-box radius 
is R = 3.5 degrees [30]. The centre of the field of view of the camera that which recorded the burst 
pointed at a = 01^01 m 56 s .49 5 = +19°20 / 38".88. The im age size is 10° x 10°. Animation available here: 
http:/ /observ.pereplet. ru/images/ GRB081130B/ grb-film. html 



The long gamma-ray burst GRB090424 

(fig USD was recorded by Swift observatory on April, 24th, 2009 14:12:09.33 UT [23J. At that time 
MASTER VWF-2 in Irkutsk operated in the alert mode (with cameras pointed apart) and GRB090424 
happened to be at the centre of the field of view of one of the cameras. It was also a typical long 
(T90 = 48sec) gamma-ray burst and UVOT detected afterglow with unfiltered magnitude of about 15.3 m 
167 seconds after the burst [24] . MASTER VWF-2 observed the grb error-box without time gaps for 1.2 
hours before, during and 1.5 hours after the gamma-ray burst with 1-second exposures. Unfortunately, 
weather conditions were far from favourable, resulting in a not so high upper limit of V\ sec < 8.0 m and 
V 60 sec < 60.0 m for a coadd of 60 images. [16J and [17J. According to both Swift/BAT [23] data and FERMI 
observatory data converted to the 15 — lbOkeV energy interval [27], 6090424 = 2.1e - 5, implying that 

7^7^090424 1 
p opt 1 



<^7T (7) 



^|r6090424 ^ 53Q 

Note that despite the rather modest upper limit the burst was very powerful and even such an upper 
limit restricts substantially the possible flux. 
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Fig. 18. Transient object found simultaneously with GRB081130B at a = 01^10 m 20 s 8 = +18°44 / 35", 14 
degrees from the centre of its last GRB081130B localisation. A detailed analysis showed the object to be a short 
flash produced by the "MOLNIYA" high-orbit satellite. 



Gamma-ray bursts recorded by FERMI observatory 

Unlike Swift, the FERMI observatory does not specialise in gamma-ray burst observations. Its onboard 
equipment allows the burst coordinates to be determined only to within several degrees or even worse. 
FERMI currently issues two to three times more alerts than Swift and is capable of recording harder 
(up to 10 MeV) radiation (Swift has an upper sensitivity limit of about 300 keV), and responds to short 
hard gamma-ray bursts, whose prompt emission has so far never been detected. Given the large error of 
the coordinate measurement, gamma-ray bursts detected by FERMI observatory may at present (and, 
very likely, also in the future) have their prompt emission observed only with wide-field cameras, such as 
MASTER- VWF. 

Until now (July, 1 2009), a total of five FERMI GRB error boxes have overlapped with the MASTER- 
VWF field of view, and two more (GRB081215 and GRB090526) proved to be outside the field of view 
after a very substantial final final refinement of their coordinates. 

GRB 08 1130 

was the first burst recorded by FERMI and synchronously observed by the MASTER VWF4 facility in 
Kislovodsk. Its coordinates and error-box have been refined repeatedly (with a characteristic scatter of 
±10 degrees), but always fell within the field of view of our cameras. Figure [T71 shows the final localisation 
of this gamma-ray burst. The error box of grb081130 has been observed without time gaps for several 
hours before, during, and for hours after the grb. 

No proper candidates could have been found inside the error box. However, given the rather large 
uncertainty of the coordinates measured by FERMI, we extended the search area for transient candidates 
to the entire frame. As a result, a very interesting object resembling a grb optical counterpart has 
been found 14 degrees from the localisation centre synchronously with grb (fig. [T8|) . However, after 
cross identification with the catalogue of artificial satellites this transient was found to coincide with 
the high-orbital satellite, which occasionally produces such flashes. Thus a rather high upper limit of 
K/r608ii30B > 12 m .O has been obtained again after the coadd of four images (X^osmos = 15 — 20s taken 
in different channels [30] ) [20] . 

When converted to the 15— 150/ceF energy interval, the gamma-ray fluence from GRB081130 is 
F srb08ii30B = LQ1 . io- 6 er#/cra 2 . We then apply the same procedures as in the case of GRB081102 

to infer FjjL 1 = 3.82 • 10~ 9 erg/cm 2 , implying a flux ratio of: g^osnsoB < ^ [33] 
090305B 

The gamma-ray burst GRB090305 (Fig. HSJ was recorded by FERMI on March, 5th, 2009 14:12:09.33 UT 
[25] . MASTER VWF-4 in Kislovodsk operated in the alert mode (with cameras apart) and observed 80 % 
part of the GRB090305 error-box at the edge of the field of view of one of cameras. Despite of its rather 
short duration (T90 = 2sec) [25], GRB090305 was nevertheless a rather long and bright gamma-ray burst. 
However, even the final position estimate has a large error-box (a = 10 /l 20 m , 5 = 68°06 / Risigma = 5.4°) 



18 



Authors Suppressed Due to Excessive Length 



Fig. 19. MASTER- VWF4 Kislovodsk GRB090305B follow-up observation. This is a 40x24 deg 2 full-frame image. 
The white and black circles show the error boxes corresponding to the 1 and 3 sigma FERMI error limits plus 
the 3-degree systematic error [25 (Risigma — 5Adeg). The black rectangle shows the IPN triangulation error-box 



Animation available here http://observ.pereplet.ru/images/GRB090305B/GRB090305B_full2.gij 



Table 2: Prompt grb observations 



GRB 



Satellite 



Site 



V 



Optical fluence 



Gamma- Gamma- 
ray fluence ray fluence 
(15-150keV) (8-1000keV) 



GRB990123 


BeppoSAX 
Konus-Wind 


ROTSE-I 


8.9 m 


4.9 • 10 


-7 




2.1 


io- 


-5 


3.0 


io- 


-4 


GRB080319B 


Swift 


Pi-of-the-sky 
TORTORA 


5.3 m 


3.7-10 


-6 




8.3 


io- 


-5 


4.4 


io- 


-4 


GRB081102 


Swift 


MASTER 
Kislovodsk 


<13.0 m 


< 2.8 • 


10" 


-8 


2.3 


10" 


-6 


3.6 


10" 


-6 


GRB090424 


Swift 


MASTER 
Irkursk 


< 9.7 m 


< 4.3- 


io- 


-8 


2.3 


io- 


-5 


5.2 


io- 


-5 


GRB0901130B 


FERMI 


MASTER 
Kislovodsk 


< 11.0 m 


< 4.1 • 


10" 


-9 


1.0 


io- 


-6 


2.0 


io- 


-6 


GRB090305B 


FERMI 


MASTER 


< 9.5 m 


< 2.3- 


io- 


-9 


2.7 


io- 


-7 


2.7 


io- 


-6 


(80% part) 




Kislovodsk 






















GRB090320B 


FERMI 


MASTER 


< 11.0 m 


< 1.5- 


io- 


-8 


8.5 


io- 


-7 


1.1 


io- 


-6 


(80% part) 




Kislovodsk 






















GRB090328B 


FERMI 


MASTER 


< 9.1 m 


< 1.4- 


io- 


-10 


2.1 


io- 


-7 


9.6 


io- 


-7 


(20% part) 




Irkursk 























and therefore no alert observations have been made. MASTER VWF-4 has observed the error-box of this 
GRB without time gaps for 3.5 hours before, during and for two hours after the with 1-second exposures. 
No optical transients close to Tqrb could be found and we obtain an upper-limit of V\ sec < 9.5 m [18]. 

F opt / 1 

rpgrb090305B ^ 1Q1 

GRB090320B 



The MASTER VWF-4 system in Kislovodsk observed synchronously 80 % of the GRB090320B (see 
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Fig. 20. MASTER- VWF4 Kislovodsk GRB090320B follow-up observation. This is a 40x24 deg 2 

full-frame image. The black circle shows the error box corresponding to the 3 sigma FERMI 

error limits plus the 3-degree systematic error (R — 12.6de(7)[26j. Animation available here 
\http://observ.perepletru/images/GRB090320B/GRB090320B_60sec.gij\ 



the GRB090424 gamma-ray burst indicates that even the measured optical limit may constitute a result 
of real importance and that and this burst deserves separate analysis. 

In this article we describe the methods of observation and analysis of images taken with very wide- 
field cameras, and their application to synchronous observations of prompt gamma-ray burst emission in 
the MASTER- VWF experiment. We finally recall that the renowned astrophysicist Bohdan Paczynski 
[4] pointed out in his last paper the exceptional astrophysical significance of robotic telescope networks, 
which, like MASTER, are composed of wide field searching cameras and larger telescopes. Now his and 
our ideas at last become reality. 

This work was supported by the Ministry of Science of the Russian Federation (state contract N 
02.740.11.0249 ). 
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Fig. 21. MASTER- VWF4 follow-up observations of GRB090328B in Kislovodsk. This is a 40x24 
deg 2 full-frame image. The black circle shows the error box corresponding to the 3 sigma FERMI 
error limits plus the 3-degree systematic error (R — 12.6deg) 21 . Animation is available here 
http://observ.perepletru/images/GRB090320B/GRB090320B_60sec7gi^ 
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